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ABSTRACT 

The SFI++ consists of ~5000 spiral galaxies which have measurements suitable for the application 
of the I-band TuUy-Fisher (TF) relation. This sample builds on the SCI and SFI samples published 
in the 1990s but includes significant amounts of new data as well as improved methods for parameter 
determination. We derive a new I-band TF relation from a subset of this sample which consists of 
807 galaxies in the fields of 31 nearby clusters and groups. This sample constitutes the largest ever 
available for the calibration of the TF template and extends the range of line-widths over which the 
template is reliably measured. Careful accounting is made of observational and sample biases such 
as incompleteness, finite cluster size, galaxy morphology and environment. We find evidence for a 
type-dependent TF slope which is shallower for early type than for late type spirals. The line-of-sight 
cluster peculiar velocity dispersion is measured for the sample of 31 clusters. This value is directly 
related to the spectrum of initial density fluctuations and thus provides an independent verification 
of the best fit WMAP cosmology and an estimate of fi^'^CTs — 0.52 ± 0.06. We also provide an 
independent measure of the TF zeropoint using 17 galaxies in the SFI-I — h sample for which Cepheid 
distances are available. In combination with the "basket of clusters" template relation these calibrator 
galaxies provide a measure of iJp = 74± 2 (random) ±6 (systematic) km s~^ Mpc~^. 
Subject headings: cosmological parameters — distance scale — galaxies: clusters: general — galaxies: 
distances and redshifts — galaxies: fundamental parameters 



1. INTRODUCTION 



The TuUy-Fisher relation (jTullv &: Fisher Ill977t ) is an 
empirical relationship between the rotational velocity 
and absolute magnitude (or luminosity) of spiral galax- 
ies. The application of this relation to measure distances 
to spiral galaxies has had an enormous impact on ex- 
tragalactic astronomy and cosmology since its discovery. 
The TuUy-Fisher (TF) relation for spirals (along with the 
Dn — <y relation for ellipticals) dramatically increased the 
volume in which redshift-independent distances could be 
measured. The TF relation and its scatter can also be 
used to constrain models of disk galaxy formation, and 
by studying the relations variation with environment, 
galaxy type and its evolution over the age of the uni- 
verse many clues into how the disks of galaxies are as- 
sembled have been found. A reliable template relation is 
an essential starting point to all this, without which no 
trust can be placed in other conclusions. Like all empir- 
ically determined scaling relations, the selection effects 
and make-up of a TF sample can introduce non-trivial 
biases into the derived template which if not taken care of 
could introduce spurious results. With improved statis- 
tics and parameter measurement both in galaxy surveys 
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and simulations the need for an accurate, unbiased tem- 
plate grows even stronger. 

The TF relation has historically played an important 
role in cosmological parameter determination, in partic- 
ular by allowing measurements of to be extended 
out of the very local universe. In fact a review of TF 
measurements of over time shows that the technique 
has done surprisingly well. The value of from TF 
studies played a significant role in the debate between 
high and low values of Hq, with cluster peculiar veloci- 
ties, the Malmquist bias and even dust extinction inter- 
vening to provide a range of measures of Hq from TF. 
The Hubble Key Project was initiated to resolve this de- 
bate and found Hq = 71 ± 2 ± 7 km s~^ Mpc~^ from 
a combination of v arious TF s urveys and Ceph e id data 
(|Sakai et al. 1 12000[ ) . Similarly iGiovanelli et al. I (|1997aD 
found Hq = 69 ± 5 km s~^ Mpc~^ from a combination 
of the I-band TF relations i n 24 clusters an d a small set 
of Cepheid cahbrators; and iTuUv fc Piercel ([2000,) found 
Hq — 77 ±8 km s~^ Mpc~^ from their sample of 5 clus- 
ters. In Section 7 below we derive Hq = 74±2±6 km s""'^ 
Mpc~^ from the sample discussed in this paper. Even in 
todays era of "precision cosmology" the TF constraints 
on Hq have similar random error to the best avail- 
able measurements combining WMAP and large scale 
structure data {e.g. Hq = 74 ± 2 km s ~^ Mpc~^ from 
a com bination of WMAP and 2dFRGS; ISanchez et al. I 
120061 ) ■ and better than WMAP alone which measures 
Hq = 73 ± 3 km s ~^ Mpc^^ (with the assumption that 
the universe is flat :ISpergel et al. |[2006l ). 

iTullv fc Fisher I (|1977D first calibrated the TF relation- 
ship in the B band using a sample consisting of 10 galax- 
ies with D < 8 Mpc and eight galaxies in the Virgo clus- 
ter. They measured M-5logh = -19.13-6.25(logVK- 
2.5), where M is the absolute magnitude, W = 2vynax is 
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the width of the Doppler broadened spectral line and 
h — HoJlOO km s^^ Mpc~^. This is equivalent to 
L (X w,^;ax- band, extinction (both internal 

to the galaxies and from our Galaxy) plays an impor- 
tant role in add i ng to the scatter of the TF relation. 
lAaronson et al. I ()1979( l argued that the relationship in 
the NIR would have a smaller scatter. They observed 
that the slope steepened in the H band, more nearly ap- 
proximating the expected L oc i^^^x from simple physical 
arguments. Much of the most succ essful early work using 
TF used the lAaronson et ai] (|198Q) calibration based on 
the distance to M31 and M33, and the assumption that 
L cx i^niax- The uncertainty in these (and other) early TF 
templates from small samples of nearby calibrator galax- 
ies was obviously dominated by small number statistics 
and since the samples did not cover a wide dynamic range 
of galaxy masses they also provided little leverage on the 
measurement of the TF slope. Several immediate at- 
tempts were made to extend the calibration, both by 
looking at cluster samples (where all galaxies could be 
assumed to be at roughly the same distance) and ob- 
taining more local calibrators. In the late 1980s it was 
argued that the Tband TF re lation could improve th e 
scatter in the relation further (|Bothun fc Mould"! 1 19871 ). 
This band, intermediate between B and H, retained most 
of the advantages gained by moving into the NIR but had 
the additional advantage of allowing CCD photometry. 

It should be obvious that using a strictly magnitude- 
limited sample to calibrate the TF relation will result 
in an underestimate of the slope, as at the low-width 
end only galaxies scattered above the line will be in- 
cluded. An early solution to this potentially strong 
bias was the adoption of so-called inverse fits in which 
the magnitude was treated as the independent variable. 
Such a fit would not suffer from incompleteness bias in 
the simple case where selection depended only on ab- 
solute magnitude. Unfortunately, realistic samples do 
not have such simple selection criteria, and in partic- 
ular internal extinction corrections, which have been 
shown to depend on galaxy luminosity, and therefore 
implicitly on the rotation wid th (Giovan elli et o/. 1119941 : 
[Masters. Giovanelli fc Havnes I (20031) make bias correc- 
tions more complicated. Arguments also exist in the 
literature that cluster samples constitute complete or 
volume-limited samples and therefore the TF relation de- 
rived from such a sample will not suffer from i ncomplete- 
ness bias. As first discussed in lSandage et ail (J 1995) this 
is not the case. An explicit or implicit magnitude limit 
in a cluster sample will have a the same impact on the 
measured TF relation as a sample with a full spread of 
distances. 

A biased TF template slope can introduce several 
subtle, qualitative biases into many of the conclusions 
which might be drawn from its application. For example 
Malmquist bias-like effects mean that the more distant 
objects in any given sample are more likely to be at the 
large-width end of the relation while the most nearby 
objects are more likely to be at the low-width end of 
the sample. If the template is biased shallow, this means 
that distant objects will preferentially have spurious neg- 
ative magnitude offsets (i.e. they appear to be artificially 
brighter than the biased template), while nearby objects 
will preferentially have spurious positive offsets. If these 
offsets were interpreted as being due to the peculiar ve- 



locities of the galaxies, this would produce a spurious 
infall region in the local universe and excess expansion 
in the more distant parts of the sample. Determinations 
of the morphological dependence of the TF relation can 
be impacted by bias on the slope too. Early type spirals 
(Sas) are preferentially brighter and have larger rotation 
widths than later type spirals. A TF template biased 
shallow might then lead to the conclusion that the zero- 
point for earlier type spirals was brighter than that for 
late type spirals. In a high redshift sample which will 
only have galaxies at the large width end of the local re- 
lation, comparison with such a biased local relation could 
easily lead to a spurious determination that the zeropoint 
of the TF relation brightens with redshift. 

In fact most high redshift stu dies of the TF re l ation use 
the relatively small samples of [ Pierce fc Tullv I (p9^ or 
iMathewson. Ford fc Buchhorn 
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as a z = com- 
B, R and I-band 
calibration is based on the TF relation of only 6 local 
spirals for which Cepheid distances were available, in 
combination with a sm all sample (~ 30 galaxies) from 
the U rsa Major cluster. IMathewson. Ford fc Buchhorn I 
()1992f ) use a sample of only 14 galaxies in the Fornax 
cluster to provide an I-band calibration. No attempt 
is made to account for sample selection biases in either 
sample. As discussed in Metevier et al. (2006), larger 
and better calib r ated local samples already exist. The 
iTullv fc Pierce I (|2000D calibration uses a significantly 
larger sample (155 galaxies in 5 clusters at I-band; 91 
galaxies at B and R and 65 at B-band), however the work 
still makes no attempt to account for selection biases 
relying on the erroneous idea that cluster samp les are 
complete. iKannappan. Fabricant fc Franx I (|2002| ) use a 
sample of 196 nearby field galaxies but just fit the inverse 
U, B and R-band relations in an attempt to account for 
bias, and make no attempt to account for the effect of pe- 
culiar velocities which also have the potential to bias the 
template slope. No fully bias corrected local template 
exists in optical wavelengths. 

The first comprehensive discussion of the effect of in- 
completeness biason the TF template relation came 
from IWillick I (fl99i ). Here it was explicitly demon- 
strated that the inverse TF relation is not bias free. A 
method was introduced which used Monte Carlo simu- 
lations to calculate the bias based on the assumption of 
a global linear TF relation, a fixed scatter and a single 
magnitude limit. Giovanelli et al. (1997b; hereafter 
G79b) extended this technique to account for a TF scat- 
ter which varied with line width, and to provide realis- 
tic accounting of sample selection criteria. They used 
this to calibrate I-band TF relation with a sample of 555 
spiral galaxies in the fields of 24 clusters. They found 
M - 51og/i = -21.01 - 7.68(logW^- 2.5), or L cx v^-^. 
They observed that the scatter in the TuUy-Fisher re- 
lation varies with the rotational velocity of the galaxy, 
from about 0.4 mag at the low width end, to less than 
0.3 mag at the high width end, resulting in distance er- 
rors for individual galaxies of 15-20%. 

Much of this paper is devoted to the derivation of a new 
template TF relation in the I-band. We follow a method 
similar to that used in G97b and use a sample of galaxies 
which builds on that cluster sample, but has a significant 
amount of new data (including the addition of several 
new clusters) making a total sample of 807 galaxies in the 
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fields of 31 clusters. This sample constitutes the largest 
ever available for the calibration of the TF relation, and 
extends the range of the calibration between Wmax = 50- 
360 km s~^ (logVF = 2.0-2.9). In combination with the 
much larger field sample (Springob et al. 2007) this work 
will finally allow for studies of the extent, non-linearity 
and morphological dependence of the TF relation at high 
statistical significance. 

This new Tband template is warranted both because 
of the addition of significant amounts of new data to the 
sample and because of changes in corrections to the raw 
data. In particular, new corrections for rotation veloci- 
ties derived from observations of both the 21cm and Ha 
lines of neutral hydrogen (HI and Ha) are expected to 
have an impact on the derived template. New simula- 
tions have been performed to study the instrument al ef- 
fects on measuring HI widths (jSpringob et al. |[2005D , and 
much greater care has been taken in combining HI and 
Ha rotation measures (Catinella, Haynes & Giovanelli 
2006). The addition of significant amount of new data 
also alters the completeness and morphological distribu- 
tion of the sample, which must be taken into account in 
rederiving the template. 

In Section 2 we discuss the observables needed for TF. 
Section 3 covers the sample selection and bias correc- 
tions. In Section 4 we discuss individual cluster TF and 
the impact of environment. Section 5 combines the in- 
dividual cluster templates into a global template. In the 
process we measure the peculiar velocity of each cluster, 
from which a cluster peculiar velocity dispersion is cal- 
culated. This quantity depends relatively simply on the 
power spectrum of initial density fluctuations and is used 
to estimate n^ '^as = 0.52±0.06. In Section 6 we describe 
the characteristics of the scatter of the TF template. Sec- 
tion 7 discuss an alternate zero-point calibration which 
uses Cepheid distances to 17 galaxies in the SFI-|--f . In 
combination with the "basket of clusters" zeropoint this 
gives a measure of Hq = 74 ± 2 ± 6 km Mpc~^. The 
final section of the paper provides a summary of the con- 
clusions. In future work including Masters et al. (2007, 
in prep.) and Spring ob et al. (2007; also see iMas ters I 
[20051: iSprineob mm this new template will be applied 
to a larger sample of spiral galaxies in the field, from 
which the peculiar velocity field, and therefore the mass 
distribution, of the local universe will be studied. 

2. TULLY-FISHER OBSERVABLES 

The sample discussed here will be referred to as the 
SFI-h-h. This sample builds on the all-sky SFI (Spi- 
ral Field I band) and SCI (Spiral Cluster I band) 
samples discussed in a serie s oilpapers in the 1990s 
(IGiovanem et al. I [l99l [l995L ilMMbl; iHavnes et al. I 
Il999allbri. b ut also includes data from the SC2 s ample 
(iDale lEggaiDale et al.l ll999D. the theses of IVogt I (|1995D 
and 'Catinella ' (|2005|) and the HI archive presented in 
[Springob et al. (200l). The entire SFI-I-I- sample wiU 
be presented in Springob et al. (2007). Subse ts of this 
data h a ve also bee n prese nted in the thes e s of [Masters I 
((2005h . ISoekkens I ((2005h and ISoringob I ([20061 ). The 
bulk of the southern hemisphe re data is presented in 
iMathewson. Ford fc Buchhorn I (1992). The fuU SFI++ 
contains ^-^5000 spiral galaxies with observations suitable 
for deriving TF distances. In this paper we consider a 
subset of that data consisting of galaxies in the vicinity 



of 31 nearby clusters. This we will refer to as the "tem- 
plate sample" of the SFI-I--I-. Most clusters in this sample 
were specifically ta rgeted for I-band observations eit her 
as£aj£of|GiovanclireraD (|1997al) . lDale et al.l (|1999f) or 
iVogt et all (2004). Fields were chosen to include at least 
one spiral galaxy for which dynamical information (ei- 
ther a HI width or optical rotation curve) was available 
at the time of observation, and the centers of the fields 
were adjusted to include as many other good TF can- 
didates (inclined, undisturbed spirals) as possible which 
might later be targeted for spectroscopy. The sample is 
therefore in no sense magnitude or angular diameter lim- 
ited, but unlike SFI/SCI does include all types of spirals. 

2.1. Optical Imaging Data 

All optical data discussed here is I-band imaging data 
taken by members of the Cornell Extragalactic group 
and their collaborators at various telesco pes. Data for 
the SFI sample of galaxies is presented in IHavnes et al. I 
(|19991j ). This was a combination of images in the North- 
ern liemis£h£re_taten^ and the sim- 
ilar [Mithewionl^rdSBuchhorn] (jl992f ) survey in the 
Southern sky from which data were carefully recalibrated 
to make a homogeneous sample. For the SFI-I--I- sample 
discussed here new data were taken in about 1900 fields 
north of (5 = -2° using the 0.9m KPNO telescope. The 
fields covered various clusters and TF candidate galaxies 
and were positioned in an attempt to optimize the num- 
ber of TF candidate galaxies observed. These imaging 
data were intended to be complete to a diameter limit 
of about 0.5 arcmin in I band, and more than double 
the number of TF candidate galaxies for which I-band 
imaging data is available. 

The Tband image s were reduced as described in 
IHavnes et all (jl999bD to provide the observed total ap- 
parent magnitude and axial ratio of the galaxy. The ob- 
served axial ratio, a/b of the galaxy comes from isophote 
fitting and is corrected for seeing as in G97b. All magni- 
tudes are co rrected Galactic extin ction using the DIRBE 
dust maps ([Schlegel et al. [[TggSD . including those mag- 
nitudes from the SFI and M FB imaging data which in 
previous publications used the lBurstein fc Heiles I (|1978l ) 
correction. Extinction due to dust internal to the 
galaxy itself is accounted for using the e xtinction correc- 
tion derived in I Giovanelli et oTI ([19941 ). A small ty pe- 
dependent k-correction as discussed in iHan I l[1992t l is 
also applied. 

2.2. Rotation Widths 

Rotation widths for the SFI-|--f sample come from a 
mixture of 21cm line (HI) global profiles and optical ro- 
tation curves (ORCs). Where they are available and of 
good quality, HI global profile widths are used in pref- 
erence to ORCs. The HI extent of a normal galaxy is 
typically about twice its optical size, so HI global pro- 
files are capable of tracing the rotation of galaxies to 
larger radii, and also have the advantage of being in- 
dependent of the assumed position angle of the galaxy. 
Approximately 60% of the galaxies in the SFI-I--I- have 
their rotational velocities measured using HI, the other 
40% use ORCs. Just less than 20% have measurements 
available from both methods. A comparison of widths 
from the different methods is presented in Catinella et 
al. (2007). 
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2.2.1. HI Spectroscopy 

The 21cm HI spectroscopy data fo r the SFI++ 
is a co mbina tion of data pu b lished in iHavnes et al. I 
(Il999af) and ISpringqb et al. I (|2005D . reprocessed by 
ISpringob et al. I ( 20051 ). For the purposes of the TF re- 
lation we are interested only in the systematic veloci- 
ties and velocit y widt hs, as discussed in Section 3.2 of 
ISpringob et al. I (|2005D . Raw widths are measured using 
an algorithm that fits low order polynomials to either side 
of the HI profile and then calculates the width at 50% of 
fp — rms (the peak flux minus the rms) on either side 
of the HI profile. These raw widths then have to be cor- 
rected for various instrumental effects, A^, the impact of 
turbulent motions At , the disk inclination, and also rela- 
tivistic effects from the observed redshift which broaden 
the observed profile by a factor {1 + z). Where the 



HI wid ths used here differ from those in IHavnes et 
(|l999a) and earlier work of this group, is in the mag- 
nitude an d trend with S/N of th e instrumental correc- 
tions, Ag. ISpringob et al. I ()2005l ) present the results of 
simulations to study the impact of spectrometer reso- 
lutions, smoothing and S/N on the measured widths, 
and their prescribed corrections have been applied to all 
the HI widths used here (including the reprocessing of 
older widths by CS). This width correction As is always 
a small positive number, which in lSpringob et al. 1 ()2005f ) 
increases linearly with log(iS'/A^) between S/N = 4-12.5. 
The previous correction for instrumental broadening de- 
creased as the S/N ratio increased. For galaxies in a 
single cluster sample where S/N roughly correlates with 
rotation width this new correction should then tend to 
shallow the TF slope slightly relative to that found with 
the old correction in G97b. The correction for turbulent 
motions is a lso done in a slightly d ifferent way from pre- 
vious work. ISpringob et al. I (|2005f ) discuss the merits of 
simple linear subtraction as opposed to a subtraction in 
quadrature. At = 6.5 km s~^ is assumed. 

2.2.2. Optical Spectroscopy 

Optical rotation curves are used to fill in areas of the 
sky or galaxies not accessible to 21cm HI observations. 
For example about half of the rotation widths in the 
southern hemisphere are from Ha/ [Nil] RCs taken at 
the 2.3m telescope at the Siding Springs Observatory. 
There are also a significant number of ORCs taken us- 
ing the long slit spectrograph on the Hale 5m telescope 
at the Palomar Observatory (see e.g. Catinella, Haynes 
& Giovanelh 2005). As d iscussed in ICatinella I (|2005f ) 
and lCatinella et ail (|2005l ) there are several algorithms 
which can be used to measure rotational velocities from 
ORCs. All of the widths used in the SFI-|— I- have been 
derived by fitting a function to the folded ORG. The 
funct ion used, the "pqlyex" model ([G iovanclli fc Havnes I 
l2002tlGatineTla1l2005HCatinella et Z , 200 5) has the form 



VpEir) = Voil - e 



-r/rpE 



)(1 + ar/rps) 



(1) 



where Vq gives the amplitude of the ORG, rpE gives the 
exponential scale of the inner rise of the ORG, and a 
gives the slope of the outer part of the ORG. This func- 
tion has no physical meaning, but provides a useful em- 
pirical fit to a wide variety of ORG shapes. Widths are 
derived from this fit by taking its value at the location 
of the optical radius, Topt which is the radius enclosing 



83% of the optical light (here derived from I-band imag- 
ing). This raw width is then corrected for inclination 
and cosmological broadening in the same way as the HI 
widths. 

As shown in iGatinella et al. I (j2007) there are system- 
atic trends of the difference between widths measured 
from ORGs and HI widths. HI widths are systematically 
larger than ORG widths presumably because HI disks 
are on average twice as large as Ha disks and ORGs 
in ge neral are flat or still rising bey ond the optical ra- 
dius (ICatinella et al. |[2006 ). Catinella et al. I pOOTil also 
show that there is a slight systematic trend of this dif- 
ference with the extent of the ORG, in the sense that 
galaxies with a larger Ha extent have larger width dif- 
ferences. This is probably because these galaxies have 
larger HI extents to match their Ha extents, so the HI 
width traces further out into the halo, while we con- 
tinue to measure the ORG width at the optical radius. 
The difference between the two measures of the rotation 
width also depends strongl y on the slope of the ro tation 
curve at the optical radius. IGatinella et al. I ()2007f ) fit for 
this trend using 873 SFI-I--I- galaxies for which both high 
quality HI widths and ORCs are available and find 

„r ,„r _ / 0.899 + 0.188w/ropt for rising ORGs 
W21/W0HC - I 1 075 - 0.013w/ropt for flat ORCs, 

(2) 

where Tmax is the maximum extent of the ORG {i.e. the 
Ha extent of the galaxy) and the definition of a flat ORG 
is one in which the gradient of the rotation curve at the 
optical radius (^opt) is less than 0.5 km s^^/arcsec; rising 
ORCs have gradients greater than this. For flat ORCs 
the correction increases the width by a fixed amount of 
~7%, almost independent of rmax/?'opt- Rising ORCs 
have corrections ranging from increases of ^^15% for 
galaxies with large Ha extents, to decreases of ~7% for 
a small number of galaxies with very small Ha extent. 
Once this correction is applied there are no systematic 
trends in the difference between HI and ORG widths, and 
the widths match within a mean scatter of ~ 5%. 

In the template sample there is a small trend of increas- 
ing Ha extent with rotation velocity such that there are 
no small width galaxies with large Ha extent, although 
large width galaxies span a range of Ha extent. The 
ORG width correction therefore makes the TF slope for 
galaxies with ORG widths slightly shallower than the 
uncorrected version. It also slightly increases the zero- 
point. The correction may also have an impact on the 
cluster peculiar velocities as derived in Section 6. Since 
the sky distribution of ORG widths vs. HI widths is quite 
uneven, some of the template clusters have virtually all 
galaxies with widths from HI measurements, while others 
are dominated by galaxies with ORG widths. A system- 
atic difference in the template zeropoint for HI and ORG 
widths could therefore be absorbed in the cluster off- 
sets (or peculiar velocities) calculated in Section 6 when 
the individual cluster samples are combined. We find 
no trends in the calculated offsets with the fraction of 
galaxies in a cluster which have ORG widths after the 
correction is applied, providing reassurance that we are 
correcting for this effect properly. 

Figure [1] shows the TF relation for galaxies in the tem- 
plate sample separated into galaxies with HI widths, and 
galaxies with widths from their ORG. The data have been 
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Fig. 1. — TF diagrams for galaxies in the in+ sample. Panel 
(a) shows galaxies with widths derived from HI global profiles; 
Panel (b) shows galaxies with widths derived from optical rotation 
curves, the widths here have been corrected using Equation |2] The 
solid line in both panels shows the bivariate fit to the total sample 
(see Section (5)1, while the dashed lines show the bivariate fits to 
the sample separated by the source of the widths. The parameters 
of the fits are given in Table |4] After the correction is applied 
to difference between the relations for the two subsamples is not 
significant. 

corrected for all biases discussed in Section 3, and galax- 
ies in each cluster sample have been shifted by the appro- 
priate amount to account for the cluster peculiar velocity. 
As is also shown in Table |4] there is only a small differ- 
ence in the TF fit to the two subsamples once the ORG 
widths have been corrected using Equation[2l The differ- 
ence is AM = -0.09 ± 0.03 -I- (0.25 ± 0.24)(log W - 2.5). 
The zeropoints differ by ~ 3ct, however we consider that 
this is not important due to the uneven sky distribution 
of the clusters. 

3. SAMPLE SELECTION AND BIAS CORRECTIONS 

3.1. Sample Selection 

Clusters are chosen to be included in the TF template if 
they have more than ~7 galaxies from the SFI-f -I- data 
set in them. The Virgo cluster is not included in the 
sample due to the complicated nature of assigning galax- 
ies to its various subgroups. Clusters too close to the 
Galactic plane are also not used, because of the uncer- 
tainty in Galactic extinction corrections. The final sam- 
ple includes 22 of the 24 clusters used by G97b to de- 
rive the template for the SFI, and an additional 9 clus- 
ters for which we now have sufhcient numbers of galaxies 
with data to include. The two missing G97b clusters are 
A2197 and A2666. Galaxies in AbeU 2197 are absorbed 
into the periphery its partner cluster A2199 instead of 
being treated separately. A2666 is a small cluster, close 
in position but at slightly lower redshift than the tem- 
plate cluster A2634; we choose not to include 9 galaxies 
associated with A2666 because of the difficulty of dis- 
entangling them from the A2634 foreground or periph- 
eral members. The assignment of galaxies to all clus- 
ters is discussed in Springob et al. (2007). We define 
two subsamples; 486 of the galaxies are considered bona 
fide cluster members, and are therefore assigned to the 
"in" sample. An additional 321 galaxies are considered 
peripheral cluster members and are added in what we 
call the "in+" sample. In deriving absolute magnitudes, 
we use the cluster redshifts for galaxies in the in sam- 
ple but the galaxy's own redshifts for those in the in+. 
Galaxies in the in+ sample generally have recessional 
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Fig. 2. — The distribution of template clusters shown projected 
onto the supergalactic plane. Positions are shown in redshift space 
(CMB frame) with no correction for peculiar velocities. Clusters 
are shown as open stars if they lie above the SOP (i.e. at +SGZ) 
and filled hexagons if they lie below it. Clusters used to define the 
"reference frame" in Section 5 are circled. The Pisces-Perseus (PP) 
region contains the nearby clusters NGC 383, NGC 507 and AbeU 
262. In the "GA" region the five template clusters are Antlia, 
Hydra, NGC 3557, Centaurus 30 and ESO 508. The concentric 
dashed circles shown radii of 3000, 6000 and 9000 km s"!. The 
Galactic plane (~ ±20°) is indicated. 

velocities very similar to that of the cluster so the choice 
between using the cluster or galaxy redshift makes little 
difference. An exception to this rule is the Cancer clus- 
ter in+ sample. Cancer is a complicated cluster with 
several distinct subgroups. Following G97b we consider 
only galaxies in the A clump of Cancer to be members 
of the in sample. Galaxies in the B, C, and D clumps 
are included in the in+ sample, and individual galaxy 
magnitudes are calculated using the recessional velocity 
of the relevant clump. 

Figure [2] shows the distribution of clusters used in the 
template in the supergalactic plane. The positions of the 
clusters are shown in redshift space (CMB frame) with 
no correction for peculiar velocity. Clusters are shown 
as open stars if they lie above the plane [i.e. at -l-SGZ) 
and filled hexagons if they lie below it. Clusters used to 
define the "reference frame" in Section [5] are circled. 

3.2. Incompleteness Bias 

As discussed in detail in G97b, the fact that all the 
galaxies in a cluster are at roughly the same distance, 
does not mean that a cluster sample is volume limited. 
Any implicit or explicit apparent magnitude limit in the 
sample means that in a given cluster we preferentially 
observe the brighter galaxies. Since the TF relation has 
an intrinsic scatter, this effect will tend to shallow the 
observed slope, by preferentially selecting galaxies which 
have brighter than expected magnitudes at the small 
width end of the relation. Similarly, it also brightens 
the zero point and decreases the observed scatter. In our 
treatment of this bias we follow G97b, and the reader is 
referred there for further explanation of the construction 
of completeness histograms for each of cluster samples 
and the Monte Carlo simulations used to derived the 
incompleteness bias from them. The completeness his- 
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tograms for each cluster are shown in Figure[3](note that 
in this Figure and elsewhere in the paper the clusters are 
listed in order of RA, and separated into Northern and 
Southern hemisphere objects). A function of the form 



is fit to these histograms; the parameters of these fits are 
shown in Table [1] and it is these functions which are used 
as inputs for the bias calculations. The final incomplete- 
ness bias corrections for the galaxies in each cluster are 
shown in Figure 31 

The incompleteness correction derived here relies on a 
priori knowledge of the intrinsic TF relation and scat- 
ter, and as such proceeds in an iterative manner. The 
magnitude of the assumed scatter has by far the biggest 
impact (when the scatter is larger the bias corrections 
are larger and vice versa), so the derived corrections 
change very little once the intrinsic scatter has been 
determined. The impact that the various assumptions 
have on the derived corrections is illustrated in Figures 
7, 8 and 9 of G97b. As well as requiring an input TF 
relation, the derivation of this bias correction must as- 
sume knowled ge of the spiral luminosity function (LF). 
As shown in Binggeli, Sandage & Tammann ' fl988), the 
spiral LF has an almost Gaussian shape. Following G97b 
we use a Schechter function with a steep faint end slope, 
a = —0.5, and M^, = —21.6. This assumption impacts 
the derived bias corrections because it determines the 
level of incompleteness in the sample. If the LF slope is 
in actual fact shallower than this the bias will be under- 
estimated (and vice versa). 

We perform several iterations of the bias correction, 
until the final TF relation differs by less than la from the 
input values. The final input TF relation is M — 5 log h = 
-20.85 - 7.81 (log - 2.5) with intrinsic scatter of e = 

0. 36-0.34(logVF-2.5). The final result is M-5\ogh = 
-20.85-7.85(log VF-2.5) and e = 0.35-0.37(logty-2.5) 
as discussed in Sections 5 and 6. 

3.3. Cluster size bias 

In deriving the template we are explicitly assuming 
that all galaxies in a given cluster are at the same dis- 
tance. This assumption varies in its accuracy based on 
the actual line of sight depth of the cluster relative to 
its distance The finite size of the cluster introduces two 
separate biases into the derived TF relation: 

1. The mean distance to a cluster. Fitting a TF relation 
under the assumption that all galaxies in a cluster are at 
the same distance, is equivalent to fitting the mean of the 
log of the distance to the cluster (because we work with 
distance modulus rather than distance). The log of the 
average is not the same as the average of the log, so this 
assumption introduces a small bias into the zeropoint of 

/3sizc = Mtrue - Mcalc = " ^ ' (4) 

where ad represents the line of sight depth of the cluster 
at a distance d. Table [T] lists the size of this bias for the 
clusters in our sample, using the assumption that they 
all have a line of sight extent of 1 Abell radius {Ra, 1-5 
Mpc/h or 150 km s^^). This is a reasonable approxima- 
tion for most of the in sample. For the in+ sample. 



we use twice this extent, increasing the small bias by a 
factor of four. 

2. Cluster size - sample incompleteness. If galaxies are 
close in magnitude to the completeness limit of the sam- 
ple they will preferentially be seen only in the foreground 
parts of the cluster {i.e. the closest parts) and therefore 
their magnitudes will be biased towards brighter values. 
This bias is illustrated in Figure 14 of G97b, and is es- 
timated assuming a Gaussian extent for the cluster, and 
using the completeness functions derived in Section 13.21 
This bias is very small (less than 0.1 mag) for most galax- 
ies in our sample. 

3.4. Edge of Catalog Bias and Homogeneous Malmquist 

Bias 

Following G97b we make no correction for either the 
edge of catalog bias or the Malmquist bias, which they 
find to be negligible in the SCI sample on which this 
sample is based. 

3.5. Morphological Dependence of the TF Relation 

Figure [5] shows TF diagrams of galaxies in the in+ 
sample separated into three groups by their morpholog- 
ical type. These data have been corrected for all biases 
discussed above and galaxies in each cluster sample have 
been shifted by an amount AM to account for the clus- 
ter's pecuhar velocity as discussed in Section O We find 
that earlier types have a shallower TF slope than later 
types as also noted in G97b. The relation for Sa galax- 
ies is shallower than that for Sc galaxies by an amount 
2.35 ±0.43 mag/dex, while the relation for Sb galaxies is 
shallower than for Scs by an amount 0.80±0.22 mag/dex. 
(see Table |4] for the full fit parameters). We need to be 
careful in interpreting these results, because the com- 
pleteness of galaxies of the different types is different in 
such a way that the incompleteness bias may be larger 
for the earlier type galaxies (as is obvious in Figure [5]). 
If true, this would artificially shallow the slope of their 
relation, perhaps accounting for the observed difference. 
In fact in the SFI sample an effort was made to be more 
complete for later type galaxies (this constraint was re- 
laxed for the additional galaxies added to the SFI-f-|-). 
Nature also conspires to make the completeness of later 
type spirals better in a TF sample since later type galax- 
ies are more likely to be strong HI or Ha emitters mak- 
ing the measurement of rotation widths easier. How- 
ever, the luminosity function of early type spirals is such 
t hat there are intrinsically fewer low l uminosity objects 
((Binggeli. Sandage &: Tammann 111988 ^ which could also 
explain the lack of small width Sa galaxies in the sample 
without the need for greater incompleteness for earlier 
type spirals. 

The observed difference in slope between morphologi- 
cal types (which is significant after a mean incomplete- 
ness bias has been subtracted) can therefore either be 
interpreted as a real physical difference between differ- 
ent types of spiral galaxies, or as an indication that the 
incompleteness bias is an underestimate for the earlier 
type spirals. Without a more extensive study of the lu- 
minosity functions for different spiral types there is no 
easy way to disentangle the two effects. The fact that 
the biggest difference is seen at the high width end of 
the relation where the bias corrections are smallest (see 
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Flc. '■'>. y\bsolutc magnitude completeness histograms for each of the in+ cluster sample. The dashed lines show the smooth fits to the 
histograms used in our Monte Carlo simulations. The template clusters are ordered by RA and separated into Northern Hemisphere (top) 
and Southern hemisphere (bottom) objects. 



Figure 5) however argues for a physical explanation. In 
either case the difference should be corrected for in the 
template. Given the small number of Sas in the sample 
(61) and the uncertainty in deriving the slope of the TF 
relation, we will correct for this effect for all 342 Sa/Sb 
spirals together (although for reference the relations are 
given for all 3 groups separately as well). We first correct 
for a plain additive offset, found by fitting all samples 
with the Sc slope. This correction is identical to that 
found by G97b (see below). After this correction we find 
the TF relation for the combined Sa/Sb sample (bivariate 



fit) to be Af-51ogft = -20.85(2)-6.92(15)(logVK-2.5), 
shallower than the Sc relation by an amount 0.95 ± 0.21 
mag/dex. We therefore apply the following corrections 
to the magnitudes: 

• SO/Sa/Sab: -0.32 - 0.9(logVK - 2.5) mag 

• Sb: -0.10 - 0.9(logVK - 2.5) mag 

• Sbc/Sc/Scd: no correction. 

If the difference in slope for the different morphologi- 
cal types of spiral galaxies cannot be explained by differ- 
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Fig. 4. — Incompleteness bias correction for galaxies in each of the template clusters shown as a function of log W. 



ing amounts of incompleteness then it is telling us some- 
thing about the physical differences between early and 
late type spirals. It is already well known that early type 
spirals arc dimmer than late types at a given rotational 
speed, what the dependence of the TF slope on morpho- 
logical type shows is that this difference gets more pro- 
nounced the more massive the galaxies are. The slope 
of the TF relation is fairly well reproduce d in numer- 
ical s imulations of disk gala.xy fo rr nation ([Mo fc Mao I 



[2O50t iNavarro fc Steinmetz I [200l . iMo fc Mao I (120001 ) 



suggest that the slope could be shallower than the fidu- 
cial L (X (observed here for Scs) if galaxies at the high- 
width end of th e relation have systematicall y less con- 
centrated halos. INavarro fc Steinmetz I (|2000) show that 
under the assumption of constant disk mass-to-light ra- 
tio, L cx v'^^^ implies that the fraction of mass in the disk. 



/disk C!C (wrot/fA)" (whcrc VA is the circular velocity of 
the dark matter halo at the virial radius). The observed 
difference in TF slope for different morphological types 
indicates therefore that early type spirals must have a 
smaller /disk at a given rotational velocity, presumably 
because they have mass in a bulge component (we find 
L cx v^-^ for Sbc/Sc/Sd, L c x f for Sbs and L cx y'^-'^ for 
Sab/Sa/SOs). FurthermoreE avarro fc Steinmetn (j2000l ) 
shows that halo concentration must increase with /disk 
to reproduce the observed TF relation; the steeper the 
slope of the relation the more quickly the concentration 
increases with /disk- Putting this together it suggests 
that as well as being intrinsically dimmer earlier type 
spirals may have systematically less concentrated halos 
than later types at a given rotational velocity. 
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TABLE 1 

Incompleteness Function Co-efficients (se e Eq uation [3t and cluster size bias 

(see Section [331 i 



Cluster Name 


Redshift 


Nin 


Nin+ 


11 

y r 


77 


(Jfll d 


in Q^iryr^ 


in-l- l3^;r,r^ 

iij. 1 /-^S!ZC 


N383 


4865* 


10 


27 


-18.98 


0.23 


0.03 


-0.002 


-0.010 


N507 


4808* 


14 


19 


-19.43 


0.83 


0.03 


-0.002 


-0.010 


A262 


4665* 


26 


49 


-18.66 


0.22 


0.03 


-0.003 


-0.010 


A397 


9594* 


14 


15 


-19.70 


0.80 


0.02 


-0.001 


-0.002 


A400 


6934* 


16 


50 


-20.08 


0.48 


0.02 


-0.001 


-0.005 


A569 


6011* 


13 


16 


-18.22 


0.05 


0.02 


-0.001 


-0.006 


A634 


7922* 


9 


- 


-18.85 


0.72 


0.02 


-0.001 


- 


Cancer 


4939* 


17 


49 


-18.96 


0.10 


0.03 


-0.002 


-0.009 


A779 


7211* 


13 


17 


-19.04 


0.45 


0.02 


-0.001 


-0.004 


A1314 


9970* 


8 


- 


-19.29 


0.17 


0.02 


-0.001 


- 


A1367 


6735* 


32 


33 


-20.47 


0.41 


0.02 


-0.001 


-0.005 


Ursa Major 


1101 


28 


34 


-16.74 


0.43 


0.14 


-0.046 


-0.186 


Coma 


7185* 


34 


43 


-20.03 


0.11 


0.02 


-0.001 


-0.004 


A2199/7 


8996* 


8 


22 


-21.11 


0.25 


0.02 


-0.001 


-0.003 


Pegasus 


3519 


19 


30 


-17.22 


0.18 


0.04 


-0.005 


-0.018 


A2634 


8895* 


18 


22 


-21.10 


0.76 


0.02 


-0.001 


-0.003 


A2806 


7867* 


10 


- 


-19.52 


0.06 


0.02 


-0.001 


- 


A2877 


6974* 


9 


- 


-19.95 


0.21 


0.02 


-0.001 


- 


A194 


5036* 


23 


31 


-17.78 


0.19 


0.03 


-0.002 


-0.009 


Eridanus 


1536 


21 


34 


-16.87 


0.33 


0.10 


-0.024 


-0.095 


Fornax 


1321 


20 


37 


-17.12 


0.41 


0.11 


-0.032 


-0.129 


A496 


9809* 


9 




-19.36 


0.44 


0.02 


-0.001 




Antlia 


3120 


17 


41 


-20.01 


0.45 


0.05 


-0.006 


-0.023 


Hydra 


4075* 


21 


31 


-19.64 


0.19 


0.04 


-0.003 


-0.014 


NGC 3557 


3120 


7 


14 


-19.71 


0.37 


0.05 


-0.005 


-0.020 


Cen 30 


3322 


23 


36 


-18.65 


0.11 


0.05 


-0.005 


-0.020 


ESQ 508 


3196 


9 


17 


-19.87 


0.35 


0.05 


-0.006 


-0.022 


A3574 


4817* 


13 


29 


-19.56 


0.17 


0.03 


-0.002 


-0.010 


Pavo II 


4444* 


8 


25 


-18.93 


0.19 


0.03 


-0.003 


-0.011 


Pavo I 


4055 


5 


16 


-19.52 


0.50 


0.04 


-0.003 


-0.014 


MDL59 


2317 


12 


25 


-17.62 


0.20 


0.07 


-0.010 


-0.042 



Cluster used to define the rest frame for peculiar velocities in Section 5.1. 



broadening and all widths corrected for inclination and 
cosmological broadening (see Section 2.2). 

1. ORC widths corrected to HI scale (Section 2.2). HI 
widths are observed to be on average larger than 
ORC widths for the same galaxy. T his correc- 
tion is derived bv lCatinella et aL~l (|2006D using 873 
SFI-I--I- galaxies which have both HI widths and 
ORCs. 

2. Morphological correction (Section 3.5). Earlier 
type spirals are observed to be on average dimmer 
than later type spirals at the same width. We also 
notice a significant shallowing of the TF slope for 
earlier type spirals. The slope correction applied is 
the same for all spirals earlier than Sc, while the 
fixed offset varies from -0.10 mag for Sb galaxies to 
-0.32 mag for spirals earlier than Sb. 

3. Incompleteness bias correction (Section 3.2). In 
any given cluster we preferentially observe galaxies 
which are scattered above the TF slope (because 
they are brighter). This has a tendency to shallow 
the TF slope, brighten the zero point and reduce 
the scatter. The correction we apply is derived 
from Monte Carlo simulations requiring an input 
TF relation, TF scatter, and assumptions about 
the spiral LF. The calculation is done several times 
until convergence is reached between the input and 
final TF. 

4. Mean distances to cluster (Section 3.3.1). This 
small fixed bias for all galaxies in a given cluster 




2 2.2 2 A S.6 2. a 

lot; W 



Fig. 5. — TF diagrams for galaxies in the in-|- sample, (a) galax- 
ies with types earlier than Sb; (b) galaxies of type Sb; (c) galaxies 
with types later than Sb. The solid line in all three panels shows 
the bivariate fit to the in-|- sample, while the dashed lines show 
the bivariate fits to the sample separated by morphological type. 
The slopes vary from -5.52 for types earlier than Sb to —7.87 for 
types later than Sb. The parameters of the fits are given in Table|4] 
Open symbols show the magnitudes corrected for all biases except 
morphological type, while the filled symbols show the magnitudes 
after the morphological type correction has been applied. 

3.6. Summary of all Bias Corrections 

Here we summarize for clarity all bias corrections 
which are applied to the data in advance of calculat- 
ing the final TF relations and scatter. These are correc- 
tions applied after the magnitudes have been corrected 
for Galactic and internal extinction and the inclinations 
have been corrected for seeing (see Section 2.1), and af- 
ter the HI widths have been corrected for instrumental 
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arises from the fact that the log of the average of 
the distance is not equal to the average of the log 
of the distance (which is what we implicitly fit for 
when using magnitudes). This bias is easy to cal- 
culate once a line of sight depth for each cluster 
has been estimated. 

5. Cluster size - incompleteness (Section 3.3.2). In 
a given cluster we preferentially observe dimmer 
galaxies only in the closest parts of the cluster - 
creating a small bias in the observed magnitudes 
which is easily corrected for using the same as- 
sumed completeness as in Section 3.2 and assump- 
tions about the size of the cluster. 

6. Edge of catalog bias and Malmquist bias. We ar- 
gue that both of these biases are negligible in the 
sample and therefore make no correction for them. 

7. Cluster peculiar velocity (Section 5). To combine 
all the individual cluster samples into a global TF 
relation, a peculiar velocity must be estimated for 
each cluster. This is done by minimizing the scatter 
in the TF relation, and the absolute value of the 
zeropoint is set by assuming that the net motion 
of a subset of the most distant clusters is equal to 
zero. 

4. INDIVIDUAL CLUSTER TF RELATIONS AND THE 
IMPACT OF ENVIRONMENT 

Here we fit the TF relation, first to each template clus- 
ter sample individually. We fit "direct" , "inverse" , "bi- 
variate" forms of the linear TF relation. The difference 
between the three linear fits is in which variable is consid- 
ered to be independent and the treatment of the errors. 
For more details of the fitting procedure see G97b. 

The individual template cluster TF relations are shown 
in Figure [5] . Galaxies in the in and in+ samples of each 
cluster are plotted as solid and open circles respectively. 
The same line showing the bivariate fit to the combined 
sample is plotted over each of the cluster sub-plots. The 
parameters of the fits to individual clusters samples are 
given in Table [2l 

4.1. Effect of Environment on the TF relation 

It is well known that environment has a strong im- 
pact on galaxy evolution. Gas rich galaxies which fall 
into clusters are thought to be stripped of their gas and 
dynamically disturbed (be it by ram pressure stripping, 
galaxy-galaxy interactions or the tidal effects of the clus- 
ter potential), therefore it is possible that the TF relation 
of galaxies in clusters may differ from that in the field. It 
should be noted here that in correcting for the observed 
morphological type dependence of the TF relation we 
may already have removed much of the impact of en- 
vironment. The observed morphology-density relations 
shows that early type spirals (which we observe to have 
a shallower TF slope than late types) are more likely to 
be found in more over-dense regions. 

We perform two simple tests to search for the effect of 
environment in our sample after correcting for morpho- 
logical type. The first is to compare the slope, zeropoint 
and scatter of individual cluster samples as a function 
of global cluster properties. The X-ray selected cluster 



sample HIFLUGCS (jReiprich fc Bohringer I l2002f ) con- 
tains 12 of the 31 clusters in the template sample. For 
these clusters no significant trends are observed of the dif- 
ference between the cluster and global slope or zeropoint 
with X-ray luminosity. There appears to be a slight ten- 
dency for clusters with higher Ljc to have smaller than 
average scatter in their TF relation, however this may 
be a result of a Malmquist bias like effect. The high Lx 
clusters are preferentially at large distances and therefore 
have a relatively small number of galaxies with TF mea- 
surements. It is observed here the TF scatter decreases 
as the number of galaxies in a cluster sample decreases, 
indicating non-Gaussian behavior in the scatter. As a 
second test for environmental effects we plot the resid- 
uals from the global TF relation as a function of the 
galaxies' projected cluster distance (Figure 7). Again no 
significant trends are observed. Both of these tests sug- 
gests that there is no reason to doubt that a global TF 
relation applies to all spirals galaxies, regardless of en- 
vironment, once a correction for morphological type has 
been applied. 

5. A GLOBAL TF TEMPLATE RELATION AND THE 
CLUSTER PECULIAR VELOCITY DISPERSION 

Here we combine the data from all 31 clusters into a 
global TF template. This is done by shifting the points 
for each individual cluster by an amount 

AAf = 51og('™^^V (5) 

V "CMB / 

taking into account the motion of each cluster in the 
CMB frame (i.e. its peculiar velocity). These offsets 
are found by assuming that all points fall on a linear TF 
relation with a global slope, and minimizing with re- 
spect the offsets. The global slope is assumed to be the 
bivariate fit to the combined sample and is found by it- 
erating until the fit converges. Additionally an inertial 
reference frame is defined by requiring that the net mo- 
tion of a subset of the most distant clusters (all those 
with cz > 4000 km s~^, excluding the Pavo Group which 
has only a small number of galaxies) is zero, i.e., 

^A^fcAAffc/^iVfc = (6) 

fc fc 

This sum is weighted by the number of galaxies in the 
cluster TF sample. The total number of galaxies in the 
reference clusters is 520 in the in+ sample, or 325 in in. 

5.1. Cluster Peculiar Velocity Dispersion 

The motions of clusters probe the peculiar velocity field 
on a scale which is well within the linear regime where 
theory and observation can be linked relatively easily. 
The large number of galaxies in each cluster also allow 
the peculiar velocities of clusters to be calculated to rel- 
atively high accuracy. Cluster peculiar velocities there- 
fore have great potential to provide constraints on cosmo- 
logical parameters. Peculiar velocities for the template 
clusters calculated from combining their TF relations are 
shown in Table [31 In that table, (as well as in Table 1) 
clusters used to define the rest frame are marked with 
a * next to their CMB velocity. Three sets of peculiar 
velocities are reported, for both the in and in+ sam- 
ples. Wpoc(l) refers to offsets when no incompleteness 
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TABLE 2 

TULLY-FlSHER RELATION FITS FOR INDIVIDUAL CLUSTER SAMPLES 



Sample 


N 




bdir 




binv 




bbi 








^b 




o 


abs 


N'^S'i in-U 
iNOOO 111^^ 


27 


-20, 


,827 


-7.623 


-20, 


858 


-9 988 


-20.844 


-8.324 





,080 


0, 


,612 


0.538 


0, 


,440 


IN OOO 111 


10 


-20, 
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bias corrections have been applied. Wpcc(2) refers to the 
offsets when the simple additive version of the morpho- 
logical correction is applied (along with all other bias 
corrections) and 'ypec(3) refers to the offsets calculated 
when the width dependent morphological correction is 
applied. The error on these offsets is taken to be the 
scatter in the individual cluster TF relation divided by 
the square root of the number of g alaxies in the cluster. 

As explained in IWatkiiisI (|1997l ). the ID peculiar ve- 
locity dispersion of clusters depends relatively simply on 
cosmological parameters, a relationship coming from its 
dependence on the initial power spectrum of fluctuations. 



In a flat universe {Qm + = 1) the relation is given by 
lOOkms"^ 



V3 



(7) 



where ag is as usual the rms fluctuation on scales of 8 
Mpc//i, Q is the mass dens ity of the unive rse, and h — 
TJo/lOO km s^i Mpc"^ In Watkins I ()1997D the function 
f{flh) is wen fit by / U.5{nh)-^-°^ + 49.4, and the 
approximation ^/f ~ 10 over the range Uh = 0.2-0.5 is 
used. The third year WMAP data release quotes fi/i^ = 
0.134 ± 0.006, /t = 0-73 ± 0.03 and erg = 0.84 ± 0.06 
(jSpergel et al. I [20061 ) . implying a value for the velocity 
dispersion of cr = 239 ± 23 km s~^. 



12 



Masters et al. 



24 
-22 
-20 
-18 
-16 



N383 / 

/ 



O 



A634 / 

/ 



-24 
22 
-20 
- 18 
-16 

-24 
-22 
-20 
-18 
-16 



- / 



N5D7 ^/ 

/ 



Cancer / 

%•/ 
■.V* 

/ 



/ Aai99/7 / Pi 



+ / 



A779 / 

. /■ 
../•• 



,!° -"- / — /• — V 

r" ° 



— / 
11^ 



A397 / A400 



/ 

-- / 

1^ i / i ii II ii 



+ 

/ 

/ 



A1314 / -- A1367 ,7 . 

/ • ; / 

/■ - / 

/ I 

' - / 

■ ■ ■ _L ill 



A569 / 

/ 

/ 



/. 



UMej / 

/ 

/; 

/»: 



/ _ _ A2634 ./ _ 2 2.5 3 

/° T ■ ^° 



2 2.5 3 



2 2.5 3 

log W 



Aaaoe / 
/• 

;/ 



-24 
-22 
-20 
-18 
-16 



Pavo 11 / 



Aaa77 J 

/ 

/ 

C 

/ 

/ 



- / 



Hydra / 

/ 



K3557 



A194 /o 



— ^. — 



/ MDL 59 / _ 



7-"° 



Eridanus / 

/ 



Cer_ 30 / 
5.^ 



Fora.x / 
./.• 



V" 

- / ° 
; / i i iiiiiii 



ESO 503 / . 

/ 



A498 / 

/ 



A3574 / 

t - 
/• 



2 2.5 3 



2 2.5 3 



2 2.5 3 

log W 



Fig. 6. — TF relations for individual clusters (not corrected for cluster peculiar velocities). Overlaid is the bivariate fit to the in sample 
from Section[5]of M — 5 log h = —20.849 — 7.845{log W — 2.5). Solid circles are in sample galaxies, open circles in+. These data have been 
corrected for the incompleteness bias and other biases discussed in Section 3. 



Figure [8] shows our measured cluster peculiar velocity 
dispersion for the in+ and in samples (here the peculiar 
velocities are calculated from the mean TF offset after 
all bias corrections to the TF data have been applied). 
The peculiar velocities of the clusters are represented as 
Gaussians of equal area, with the width being set by the 
error on the measurement. An arbitrarilly scaled sum of 
all these Gaussians is shown as the heavy solid line. A 
Gaussian fitted with zero mean to the entire in+ sample 
gives a = 365 ± 34 km s~^, (or a = 440 ± 40 km s~^ in 
the in sample). Here the errors are calculated using a 
Monte Carlo simulation in which 1000 realizations of the 
cluster peculiar velocities (drawn from the normal distri- 
bution shown individually for each cluster in Figure [8]) is 
made. This measurement is biased larger than the intrin- 



sic cluster peculiar velocity dispersion because of broad- 
ening due to the error on the individual cluster peculiar 
velocity measurements. A simplistic correction for this 
broadening, (subtracting the mean peculiar velocity error 
in quadrature) gives an intrinsic cosmic cluster velocity 
dispersion of cr = 287±34 km s~^ and 365±40 km s~^ for 
the in+ and in samples respectively. Correcting for 
error broadening slightly more carefully using another 
Monte Carlo simulation gives a = 298 ± 34km s^^for 
in+ and a = 368 ± 40km s^^for in. 

Our measurements of the cosmic cluster peculiar ve- 
locity di spersion are entire l y con sistent with those calcu- 
lated bv lGiovanelh et al. I (|l998l ) of 270±54 km s"! and 
277±63 km for their in+ and in samples respectively. 
These values are also in agreement with other measures 



SFI++ TF Template 
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Distance from cluster center (Mpc/h) 



Fig. 7. — Residual from the global TF relation (bivariate fit to 
the in+ sample) for all galaxies in the in+ sample as a function 
of their projected distance from their cluster center. 

using the SCI data ijBahcall fc Oh lll99ftlWaAkhIs1[l99l . 

which aU favor aid ^ 250 km once values have been 
corrected for error broadening. 

This measurement of the cosmic cluster peculiar ve- 
locity dispersion is also consistent with the WMAP es- 
timate, providing yet another independent verification 
of the best fit cosmological model. We can use the 
approximation ^/J ~ 10 in Equation 7 to calculate a 
value n'^-^ag = 0.52 ± 0.06 from this TF measurement 
of the cluster peculiar velocity dispersion alone (where 
the error here only represents the measurement error on 
the cluster peculiar velocity dispersion and does not ac- 
count for any bias introduced by the assumption). Al- 
ternatively we can use WMAP information on ilh^ — 
0.134 ± 0.006 (coming from the relative amplitude of the 
power spectrum peaks) along with our measurement of 
h = 0.74 ± 0.02 (see Section 8) from a combination of 
TF and Cepheid data to calculate ^T/) = 11.34 ± 0.59, 
implying n°-^as = 0.46 ± 0.06. 

These cluster peculiar velocities also provide tanta- 
lizing hints about the local flow field. It is interest- 
ing to note that the two clusters, Hydra and A3574 
at cz > 4000 km s~^ in the "Great Attractor" region 
both have slightly negative peculiar velocities (A3574 has 
zero peculiar velocity to within la), while the clusters at 
cz < 3500 km s^^ in that direction have positive pe- 
culiar velocities. This is presumably the signature of in- 
fall onto the "Great Attractor" and could be interpreted 
as a detection of backside infall into that region. Infall 
onto the Pisces-Perseus supercluster is also hinted at. 
Abell 262 and NGC 383, both at D - 65 Mpc//i75, have 
very marginal negative peculiar velocities in the CMB 
frame (N383 has Wpcc — 0km s~^ to well within la) and 
are presumably falling into the supercluster from behind, 
while NGC 507 at -D ~ 60 Mpc//i75 in the same direction 
has a positive peculiar velocity indicating front-side in- 
fall. Also of note, the pair of clusters A2806 and A2877 
have opposite peculiar velocities implying very similar 
real space distance, and presumably showing their mo- 
tion towards each other and indicating a future cluster 
merger. In the direction of the constellation Pavo, Pavo 
I (or the Pavo group) at D ~50 Mpc appears to be in the 
foreground of Pavo II (or the Pavo cluster) at D ~ 60 
Mpc and moving towards it. 

5.2. The Global TF Template 




Fig. 8. — Line of sight peculiar velocities in the cluster reference 
frame all clusters. Peculiar velocities are represented as equal area 
Gaussians, with the width being set by the error on the measure- 
ment. Panel (a) shows the results using the in-j- sample for each 
cluster; Panel (b) shows the in sample. The heavy solid line in 
both panels is an arbitrarily scaled sum of all the Gaussians repre- 
senting the individual cluster peculiar velocity measurements. The 
dashed and dotted lines are fits to the sum, the dotted line having 
the mean fixed to zero. When the mean is allowed to vary it is not 
significantly different from zero. 

The global template TF relation is shown in Figure [9] 
for both the in and in+ samples. We also show the in+ 
sample plotted with error bars. A summary of linear fits 
(direct, inverse and bivariate) to various subsets of the 
full sample is shown in Table ID Quadratic fits do not 
improve the scatter around the relationship significantly 
and therefore are not reported. The first two lines in 
Table 4 give the TF fits for the in and in+ samples 
when no incompleteness bias (see Section 3.2) is applied. 
The direct and bivariate fits in these lines ought to have 
shallower slopes and brighter zeropoints than the bias 
corrected relation, which is observed. Next we show TF 
fits for incompleteness bias corrected magnitudes. The 
first set in Table 4 show the fits for different types of spi- 
ral galaxies. These fits are discussed in Section 3.5 and 
used to argue for a morphological type bias correction. 
Next we show fits using a plain additive morphological 
correction (as used in G97b, and discussed in Section 
3.5) to account for the fact that at a given rotational 
width earlier type spirals are dimmer than later types. 
We show fits for in and in+, and also for in+ divided 
into Scs and types earlier than Sc. Here it is clear that 
the slope for the earlier type spirals is shallower than the 
Sc slope (as previously discussed in Section 3.5), suggest- 
ing that a type dependent morphological is needed. This 
correction is applied to the magnitudes used in the fits 
shown in the last part of Table 4. It should be noted here 
that the addition of this width dependent morphological 
correction is in effect causing us to fit two TF relation 
to each of these sets of galaxies, one for the Scs and a 
second for galaxies earlier than Sc. Here we show fits for 
the in and in+ samples, but also the in+ sample with 
a 2.5ct clip applied to reduce the impact of outhers; for 
the in+ sample with inclination cuts applied since less 
inclined galaxies have larger errors on inclination correc- 
tions to the widths, and the most inclined galaxies can 
have large errors introduced from the impact of internal 
extinction. Finally we show the fit divided into source 
of rotation width (as discussed in Section 2.2.2), to show 
that there are no sig nificant differenc e s onc e the ORG 
correction derived in iGatinella et al. I ()2006l ) is applied 
to correct ORG widths to a HI width scale. 

The zeropoint derived for the TF relation from the 
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Fig. 9. — Global TF relation for (a) in sample, (b) in+ sample, 
(c) in+ sample plotted with measurement errors. The direct and 
bivariate fits to the respective samples are overlaid as dashed and 
solid lines respectively, while the G97b bivariate fit to their in+ 
sample is plotted as the dotted line. 

"basket of clusters" method could have a systematic off- 
set from the real value if there is a net motion of the 
cluster reference frame with respect to the CMB. We 
can estimate this offset under the assumption that the 
peculiar velocities of clusters are characterized by an 
rms ID velocity dispersion cr„. The mean velocity of 
N randomly chosen clusters with uncorrelated veloci- 
ties is then (vpcc) ~ cr«-/V~^/^. In practice all 21 of 
our reference frame clusters do not constitute uncorre- 
lated measurements of the peculiar velocity field, which 
is known to show coherent structures over large scales. 
A conservative estimate suggests we have ~ 12 indepen- 
dent measurements of the velocity field. We measured 
(Tt, — 298 ± 34 km s^^ in the preceding section, imply- 
ing (wpcc) ~ 86 km s~^ for the set of reference clusters 
(assuming there are no coherent flows on the scale of the 
entire sample volume). The mean redshift of these clus- 
ters is {cz) — 6750km s^^. The likely magnitude offset, 
\AM\ - 2.17(upec)(c^)-\ is therefore |AAf| ~ 0.03 mag. 
This value should be combined in quadrature with the 
statistical error on the zeropoint of ~ 0.015 mag for the 
in+ sample (0.02 mag for in) to provide the final zero- 
point accuracy of 0.03 mag for the in+sample (or 0.04 
mag for in). 

5.3. Comparison with SCI TF Relation 

As discussed above, the sa mple here is based on the 
cluster sample SCI presented in lGiovanelli et aC\ (|1997aD 
to which a global TF relation is fit in G97b, but includes 
a significant amount of additional data, including 9 ex- 
tra clusters (for which previously there was not sufhcient 
numbers of galaxies with TF measurements to include) 
totaling ^ 300 new galaxies (or 30% of the sample). The 
methods we have used to correct for biases, and to fit 
the TF relation are very similar to those discussed in 
G97b, so a direct comparison of the global TF relation 
should be made. The template sample here differs from 
the SCI sample mostly by the addition of small diameter 
and faint apparent magnitude galaxies, so for galaxies 



in a cluster all at basically the same distance this will 
preferentially add galaxies at the low mass end of the 
TF relation, adding extra leverage to the measurement 
of the slope, and increasing the overall scatter. As dis- 
cussed earlier, the new instrumental corrections applied 
to the HI widths (which make up about 60% of the veloc- 
ity widths used) should also shallow the TF slope slightly 
relative to what would be found with the old corrections, 
and we also apply a correction to put ORC widths onto 
the same scale as the HI widths, which slightly shal- 
lows the slope and dims the zeropoint for these galaxies. 
There is also a significant difference in the distribution of 
different morphological types in the samples. SCI was de- 
signed to be heavily dominated by Sc galaxies, with 63% 
of the template sample being Sc or later. In this tem- 
plate sample 58% of the galaxies are Sc or later, meaning 
that more than 50% of the new galaxies are earlier type 
spirals. Here we apply a width dependent morphologi- 
cal type correction to the data since with the improved 
statistics the difference in slope between morphological 
types appears significant. With this correction our final 
TF relation fit to the in+ sample is 

M-5\ogh = -20.85(2) - 7.85(10) (log - 2.5), (8) 
using the plain additive correction for morphological 
type as in G97b we find (for the same galaxies) M — 
5\ogh = -20.82(2) - 7.43(ll)(log W - 2.5). The values 
of these slopes bracket the G97b result of M — 5 log h — 
-21.01(2) -7.68(10)(logVF- 2.5), which is what should 
be expected for a sample which has more early type spi- 
rals than the SCI sample. The zeropoints of the G97b 
TF relation and our favored fit differ by an amount 
AM — 0.16±0.03, i.e. they are different at the 5a level, 
such that the G97b value is 0.16 ± 0.03 mag brighter. 
Part of this difference can be ac counted for by the use 
of DIRBE (jSchlegel et al. Ill998f ) corrections for Galac- 
tic ex tinction, which are larger than those used in G97b 
(from iBurstein fc Heiles Ill978f ) by a mean value of 0.05 
mag over lines of sight towards the template clusters, but 
this still leaves a difference of AM — 0.11 ± 0.03 mag. 
This can probably be explained by the use of the ORC 
correction from Section 2.2.2 (which was not available 
for ORC widths in the G97b template). This correction 
should dim the magnitude of the zeropoint of the whole 
sample slightly (because it increases the widths of galax- 
ies with ORCs slightly while they remain at the same 
magnitude) . Without applying it we find both a slightly 
larger cluster velocity dispersion (as discussed in Section 
2.2 part of the difference between the ORC and HI widths 
can be absorbed into the cluster peculiar velocity disper- 
sion because most clusters are biased towards the use of 
one method for measuring widths) and a zeropoint of ~ 
-20.9 mag. The difference between this zeropoint and 
the G97b value (after adjusting for the DIRBE dust cor- 
rection) is ^ 0.1 ± 0.03 mag, now only a 3a difference, 
although it should be noted here that since the samples 
are not independent, the significance of this agreement 
is not the same as for totally independent studies. 

6. THE SCATTER IN THE TF RELATION 

The intrinsic scatter in the TF relation carries as much 
(if not more) information for models of galaxy formation 
as the slope and intercept, and a proper understanding 
of the scatter is also important in deriving errors on dis- 
tances (or peculiar velocities) from TF. As discussed in 
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TABLE 3 

Fit off^et'< for ixDn'inrAL clt'ster samfles 



Sample 


"CMB 


fpec(l) (K 


:m S ^) 


i'pec(2) (km s 1) 


Vpec{3) (km s ^) 




(km s ) 


in 


in-h 




in+ 




in-l- 


N383 


4865* 


-24± 227 


-115+ 254 


-5+ 230 


-58+ 244 


-25+ 226 


-76+ 239 


N507 


4808* 


393+ 243 


236+ 199 


351+ 235 


208+ 191 


413+ 231 


253+ 190 


A262 


4665* 


_Q2S+ 21 2 


-151+ 157 


-248+ 225 


-67+ 157 


-240+ 216 


-82+ 153 


A397 


9594* 


629± 698 


443+ 687 


512+ 745 


313+ 727 


622+ 723 


423+ 705 


A400 


6934* 


196+ 217 


267+ 141 


137+ 208 


225+ 141 


150+ 210 


229+ 141 


A569 


6011* 


-438+ 227 


-436+ 249 


-155+ 200 


-193+ 228 


-157+ 195 


-215+ 229 


A634 


7922* 


-773+ 459 


-831+ 448 


-587+ 487 


-664+ 479 


-577+ 475 


-663+ 462 


Cfiiicer 


4939* 


23+ 264 


8+ 146 


203+ 256 


95+ 142 


153+ 258 


77+ 144 


A779 


7211* 


-599+ 302 


-628+ 255 


-616+ 317 


-603+ 264 


-582+ 323 


-600+ 267 


A1314 


9970* 


-107+ 514 


-189+ 520 


332+ 522 


230+ 529 


84+ 544 


-33+ 556 


A1367 


6735* 


192+ 212 


187+ 206 


32+ 221 


27+ 213 


71+ 216 


66+ 209 


Ursfi Kdcijor 


1101 


-392+ 53 


-378+ 54 


-352+ 53 


-338+ 54 


-348+ 51 


-333+ 52 




7185* 


619+ 165 


RQQ-I- 155 

Kj Hf tj -L- A-^\J 


323+ 142 


358+ 143 


361+ 142 


398+ 142 


A2199/7 


8996* 


-82+ 344 


441+ 317 


-407+ 378 


145+ 334 


-146+ 350 


251+ 316 




3519 


-317+ 135 


-450+ 131 


-164+ 134 


-314+ 131 


-131+ 139 


-290+ 130 


A2634 


8895* 


820+ 298 


540+ 269 


652+ 284 


257+ 280 


660+ 295 


318+ 277 




7867* 


Q4n-I- 273 


Qoi;-!- 270 


395+ 312 


365+ 309 


266+ 280 


236+ 279 


A2877 


6974* 


-309+ 379 


-339+ 390 


-361+ 462 


-407+ 481 


-596+ 479 


-647+ 502 


A194 


5036* 


-295+ 281 


-471+ 236 


-141+ 258 


-352+ 229 


-168+ 256 


-348+ 220 


Eridanus 


1536 


-378+ 113 


-378+ 80 


-307+ 108 


-316+ 78 


-306+ 109 


-310+ 78 


Fornax 


1321 


-226+ 72 


-260+ 59 


-167+ 74 


-202+ 57 


-185+ 72 


-209+ 56 


A496 


9809* 


-84+ 497 


-141+ 503 


72+ 501 


-6+ 512 


15+ 479 


-71+ 490 


Antlia 


3120 


290+ 105 


203+ 84 


336+ 99 


205+ 89 


292+ 102 


173+ 89 


Hydra 


4075* 


-488+ 157 


-284+ 127 


-416+ 182 


-221+ 137 


-422+ 169 


-246+ 129 


NGC3557 


3318 


227+ 167 


229+ 199 


308+ 153 


191+ 210 


281+ 160 


217+ 202 


Ccii30 


3322 


216+ 135 


278+ 100 


249+ 126 


331+ 96 




324+ 94 


ESQ 508 


3196 


411+ 129 


232+ 102 


396+ 113 


254+ 96 


382+ 122 


230+ 99 


A3574 


4817* 


-190+ 209 


-27+ 128 


-191+ 214 


-20+ 134 


-199+ 210 


-26+ 132 


Pavo II 


4444* 


302+ 179 


11+ 164 


314+ 162 


10+ 160 


304+ 147 


34+ 155 


Pavo I 


4055 


393+ 202 


173+ 164 


489+ 186 


191+ 163 


473+ 191 


212+ 159 


MDL59 


2317 


-535+ 135 


-347+ 128 


-411+ 137 


-303+ 113 


-410+ 134 


-277+ 117 


* Cluster used to define the rest frame for peculiar 


velocities in 


Section 5.1. 
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Linear TF parameters for global samples. 






Sample 


N 








O-hi hi 


Co 


^6 <^bi f^abs 



No incompleteness correction 



in 

in-l- 



486 
807 



-20.906 

-20.915 



-7.287 
-7.190 



-20.883 
-20.890 



-8.294 

-8.341 



-20.892 

-20.901 



-7.565 

-7.442 



0.019 
0.015 



0.134 0.384 0.301 
0.105 0.411 0.321 



Incompleteness correction as described in Section 3.2 



No morphological correction 



in-|-, Sa 


61 


-20.826 


-5.199 


-20.945 


-7.819 


-20.774 


-5.516 


0.073 


0.403 


0.399 


0.317 


in-l-, Sb 


281 


-20.742 


-6.805 


-20.847 


-8.506 


-20.725 


-7.072 


0.025 


0.172 


0.379 


0.300 


in-|-, Sc 


465 


-20.855 


-7.685 


-20.956 


-8.726 


-20.853 


-7.870 


0.020 


0.145 


0.417 


0.328 


Plain additive morphological correction 


in 


486 


-20.820 


-7.318 


-20.844 


-8.031 


-20.808 


-7.561 


0.019 


0.134 


0.388 


0.303 


in-l- 


807 


-20.833 


-7.192 


-20.853 


-8.087 


-20.823 


-7.425 


0.015 


0.106 


0.413 


0.324 


in-t-. Sa/Sb 


342 


-20.873 


-6.606 


-20.851 


-8.445 


-20.845 


-6.919 


0.023 


0.154 


0.392 


0.310 


in-|-, Sc 


465 


-20.855 


-7.685 


-20.956 


-8.726 


-20.853 


-7.870 


0.020 


0.145 


0.417 


0.328 


Width dependent morphological correction 


in 


486 


-20.841 


-7.765 


-20.866 


-8.409 


-20.829 


-8.006 


0.020 


0.137 


0.382 


0.299 


in-l- 


807 


-20.859 


-7.625 


-20.881 


-8.435 


-20.849 


-7.845 


0.015 


0.103 


0.407 


0.320 


in-|- 2.5iT clip 


794 


-20.857 


-7.715 


-20.898 


-8.455 


-20.845 


-7.920 


0.015 


0.098 


0.384 


0.307 


in-l- i > 45° 


745 


-20.855 


-7.710 


-20.902 


-8.485 


-20.847 


-7.868 


0.015 


0.102 


0.401 


0.315 


in-l- 60°< j < 80° 


411 


-20.876 


-7.711 


-20.922 


-8.476 


-20.873 


-7.842 


0.020 


0.136 


0.393 


0.312 


in-|-, HI widths only 


577 


-20.886 


-7.578 


-20.920 


-8.481 


-20.871 


-7.773 


0.017 


0.111 


0.412 


0.323 


in-|-, ORG widths only 


230 


-20.793 


-7.654 


-20.840 


-8.364 


-20.777 


-8.026 


0.030 


0.218 


0.371 


0.301 
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Section [3?2l the observed scatter in the TF relation of any 
sample with an implied or actual magnitude limit will be 
an underestimate of the true amount. Magnitude incom- 
pleteness in the sample means that galaxies which are 
intrinsically dimmer than the relation are preferentially 
removed when they are near the magnitude limit. This 
has the effect of shallowing the observed slope, brighten- 
ing the observed zeropoint and decreasing the observed 
scatter (especially at the small width end of the rela- 
tion). We model the bias, by making 1000 realizations 
of the global cluster sample, using the completeness for 
individual cluster samples as calculated in Section 13.21 
The result for all galaxies in the in+ sample is shown 
in Figure [TOl which also shows the average bias in bins 
of 20 galaxies, and a parametric fit which we use to cor- 
rect the measured scatter. There is a wide range in the 
amount of bias at small log W because of the wide range 
of completeness at low absolute magnitudes found in the 
individual cluster in+ samples. Above logVF = 2.8 we 
assume the bias to be negligible, and the parametric fit is 
fixed at that point. The bias correction for logM^ < 2.8 
is 

'^""'"'"''"'^ = 1.0-0.0044(logVr-2.8)-0.367(logW-2.8)2. 

ftruc 

(9) 

The measured and bias corrected scatter (standard de- 
viation) is shown by the open and solid square points in 
Figure [TTl A linear fit to the bias corrected points gives 
e^bs = 0.41 - G.44(logM^ - 2.5) (shown by the straight 
solid line in the figure). This total scatter is somewhat 
larger than eobs = 0.32 - 0.325(logW^ - 2.5) found by 
G97b for the SFI template sample, but has a similar de- 
pendence on log The source of the additional scatter 
must be related to the difference in the SFI and SFI-f + 
samples in terms of both the distribution of galaxy prop- 
erties and changes in the way measured values are cal- 
culated. The template sample of SFI-I--I- adds to the 
SFI mostly small diameter galaxies, and also more early 
type spirals than in SFI which was designed to be dom- 
inated by Sc galaxies. These small galaxies are harder 
to assign proper morphological types (adding scatter in 
the morphological correction) , and harder to measure in- 
clinations for, adding a significant amount of scatter to 
the widths. The SFI-f-t- template sample should there- 
fore both have higher measurement errors on average and 
possibly also higher intrinsic scatter. We use the entire 
SFI-f-|- sample (both template and field galaxies) to es- 
timate the total measurement error at a given rotation 
width for SFI-I--I- galaxies. These values are shown by 
the dotted lines in Figure [TT] The lower line at cr ~ 0.1 
shows the measurement error on the total magnitudes 
which at all widths is the least important source of error, 
and is similar to the value for SFI galaxies. The middle 
line shows a value of 7.85aw which is the error on the 
measured rotation widths expressed in magnitudes. This 
quantity, particularly at small widths appears larger than 
that found for SFI galaxies. Adding these two contribu- 
tions in quadrature approximates the total measurement 
error (this neglects the covariances between the two val- 
ues which arise from the inclination corrections) which 
as expected is slightly larger than the total measurement 
error found for SFI. 

The dashed lines in Figure [TT] show the total measure- 
ment error for SFI-I-+ galaxies with fixed scatters of 0.2, 





.' 1 1 1 ' ' ' Y 1 ' . 










■• . • JTh'T \--- ■ ■! 

J ■ ¥ 1 . ■.tw I- i~i 1 • ■ ■ 


- 




. ^ 

■ 1 ■ 1 I,- "1 I'l '111-.. 




J 






1,; ..^ : 




, , ■ , 1 , , , 1 , , , 1 , - 



a 2,2 2,4 2,6 2,8 



log W 

Fig. 10. — Our calculation of the underestimate of the scatter in 
the TF relation due to the incompleteness bias for each in+ clus- 
ter. The bias for each galaxy was estimated from 1000 realizations 
of the individual in+ clust er sa mples with incompleteness charac- 
teristics derived in Section 13.21 The large squares show the mean 
bias in bins of ~20 galaxies. The solid line is a parametric fit to 
the bias which we use to correct the measured scatter. 

0.25, 0.3, 0.35 and 0.4 mag added in quadrature. As is 
obvious this does not provide a good fit to the measured 
scatter. At the low width end a much larger intrinsic 
scatter is needed than at the high width end even af- 
ter the larger measurement errors have been accounted 
for. We therefore fit for an intrinsic scatter which has a 
width dependence and find eint = 0.35— 0.37(logW^— 2.5). 
This fit is shown by the dot-dashed line in Figure [TTl 
and provides a better fit to the observed error than 
the straight line fit. This value for the intrinsic scat- 
ter is again larger than that found in G97b, who measure 
eint = 0.26 -0.28(logW^- 2.5) for the SFI cluster sample. 
Even accounting for the increased measurement errors in 
SFI-f-l- there still appears to be a larger intrinsic scatter, 
presumably related to the larger variety of morphological 
types and perhaps the addition of higher redshift clusters 
into the template sample. 

7. A CEPHEID CALIBRATION OF THE TF ZEROPOINT 
AND Ho. 

In this section we use the classic "distance ladder" ap- 
proach to provide an independent calibration of the TF 
zeropoint. The "basket of clusters" method used above 
has several advantages over this classic approach, the pri- 
mary being that it is independent of the Hubble constant. 
Even in today's era of "precision" cosmology, the error on 
the Hubble constant will significantly increase the error 
on the TF peculiar velocities measured using a template 
calibrated with the distance ladder method. Peculiar ve- 
locities from the "basket of clusters" method are indepen- 
dent of Hq. a further disadvantage of the distance ladder 
method is the still small numbers of calibrators. In the 
entire SFI-|--f sample of ~ 5000 galaxies, only 17 have 
distances measured from their Cepheid variables (data 
shown in Table These galaxies tend to lie at the high 
width end of the relation which could introduce a bias on 
the zeropoint. They also tend to have relatively low incli- 
nations and thus the error on the determination of their 
actual rotation speeds is large. However, it is possible 
that there is a net motion of the reference cluster sam- 
ple (see Section [5]) in the CMB frame which would bias 
the zeropoint measured with that method. The refer- 
ence sample includes clusters with recessional velocities 
between 4000-10000 km/s, with a mean value of 6750 
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Fig. 11. — Total observed scatter in the TF relation fit to the 
in+ sample. The open squares show the observed standard de- 
viation. The filled squares have been corrected for incomplete- 
ness bias as shown the parametric fit in Figure 1101 The solid 
line shows a linear fit to this bias corrected standard deviation 
of Eobs = 0.41 - 0.44(logVy - 2.5). Also plotted is the error bud- 
get for all galaxies in the SFI-I — h dataset as a function of rotation 
width. The dotted lines show errors associated with the photom- 
etry and width measurement (multiplied by the TF slope to be 
expressed in magnitudes). The sum in quadrature of these two 
approximated the total measurement error. The dashed lines show 
a sum in quadrature of this total error with fixed scatters of 0.2, 
0.25, 0.3, 0.35 and 0.4 mag respectively. The dot-dashed line shows 
a sum in quadrature of the total scatter with a width dependent 
intrinsic scatter of e^^t = 0.35 - 0.37(log W - 2.5). 

km s""'^. These clusters are distributed fairly evenly over 
the sky, but a small anisotropy in the distribution com- 
bined with the possibility of relatively large bulk flow 
towards the Shapley supercluster at cz ~ 12000 km s^^, 
could introduce a bias. 

It therefore seems prudent to check the "basket of clus- 
ters" calibration with an independent measure of the ze- 
ropoint from the distance ladder approach. Since there 
are no galaxies in the SFI-|--f sample with W < 200 
km s~^ that have measured Cepheid distances, it is not 
appropriate to derive a slope from the calibrator sample, 
so we will use the slope derived in Section [5] and fit for 
the zeropoint alone. Under the assumption that the ze- 
ropoint calculated from the "basket of clusters" method 
does provide an unbiased measurement, this comparison 
alternatively provides an estimate of Hq. 

The bulk of the Cepheid distances for SFI-I--I- galaxies 
come from the HST Key Project (Freedman et al. 2001; 
FOl). Two additional distances come from Leonard et 
al. (2003; NGC 1637) and Newman et al. (1999; NGC 
4603, this distance has been adjusted from that reported 
in the paper to be consistent with the newer Cepheid cal- 
ibrations). We give in Table [5] the extinction corrected 
distance moduli reported in those papers. Cepheid dis- 
tances are known to ha ve a dependence on the metallic- 
ity of the host gal axy (jSakai et al. Il2004 ). We use the 
ISakai et al. I (|2004[ ) correction of (m — M)cor — {m — 

M)-(0.24±0.05)([O/H]gai-[O/H]LMc) with [0/H]lmc 
= 8.50 dex, propagating the uncertainty in this correc- 
tion through to the final Cepheid distance modulus un- 
certainty. The apparent magnitudes listed in Table [5] 
have been corrected for Galactic and internal extinction 
and the width-dependent morphological type correction 
from Section 3.5. All but two of these widths come from 
HI measurements. The two ORC widths for galaxies in 
the Fornax cluster are noted in Table [5] and small cor- 
rections from Equation [2] have been applied (both of the 



ORCs are fiat so the corrections are very small). Panel 
(a) of Figure [12] shows the galaxies with Cepheid dis- 
tances on the TuUy-Fisher relation, while their residuals 
from the bivariate fit to the in+ sample are shown in 
panel (b). While the TF relation from Section [5] is de- 
rived independent of the value of Hq, a value must be 
used to compare with the Cepheid calibrated zeropoint. 
In Figure fT2l the best fit relat ion has been adjuste d for 
iJo 74 ± 2 km s'^ Mpc"! (jSanchez et al. II2006I from 
a combination of WMAP and 2dFRGS); the upper and 
lower limits on Hq are shown by the two solid lines in 
both panels. The errors on the residuals for the individ- 
ual galaxies show the sum in quadrature of the measure- 
ment error on the apparent magnitudes, the error on the 
Cepheid distance, and the error on the measured width 
multiplied by the slope of the TF relation. The error on 
the residuals for all these galaxies is dominated by the 
measurement error on the galaxy widths. 

The error weighted mean of the residuals for these 
galaxies gives Mcophoid - Mciustors = -0.01 ±0.05 ±0.06. 
Here the first error quoted is the statistical error on the 
weighted mean, the second takes into account the 3% er- 
ror on Hq. This fit is shown by the dotted line in both 
panels of Figure 1121 It is indistinguishable from zero. 
The Cepheid calibrated TF relation (using the slope from 
the bivariate fit to the in+ sample) is therefore 

M = -20.485 ±0.05 ±0.06- (7.85 ±0.1)(logiy- 2.5), 

(10) 

(where again the first error is statistical and the second 
due to the uncertainty on Hq). This is to be compared 
to 

M-51og/i74 -20.474±0.03-(7.85±0.1)(logM^-2.5) 

(11) 

from the basket of clusters method. Here we note again 
that this template is independent of Hq, but has been 
adjusted to Hq — 74 km s"^ Mpc^^ for ease of compar- 
ison. It is certainly reassuring that the two methods for 
calibrating the zeropoint agree so well. This agreement 
provides a strong limit on how much the cluster reference 
frame differs from the CMB rest frame. 

Under the assumption that the "basket of clusters" 
method provides an unbiased measurement of the TF ze- 
ropoint, we can use the combination of that TF relation 
and the independent Cepheid calibration of the zeropoint 
to provide an estimate of Hq. Equating the zcropoints 
from the two methods is equivalent to using the Cepheid 
calibrated TF relation to measure redshift-independent 
distances for the 31 clusters in our sample and taking the 
mean of the values of Hq calculated from v/D for each 
cluster. Since we account for the cluster peculiar ve- 
locities when combining their TF relations (see Section 
5) here v is already corrected for the peculiar motion 
of each cluster. By equating the two zeropoint we derive 
Hq = 74±2 (random) km Mpc~^. This measurement 
has two possible sources of systematic error. The first is 
that the "basket of clusters" template may be measure a 
biased value for the TF zeropoint if there is a net motion 
of the reference frame. This point is discussed in more 
detail in Section 5.2, where we estimate at most a system- 
atic offset of 0.03 mag which is included in the cluster TF 
zeropoint uncertainty (and therefore in the error already 
quoted for Hq). A larger source of possible systematic 
error comes from the Cepheid distance scale. The er- 
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TABLE 5 

Galaxies in SFI++ with Cepheid Distances 



Galaxy 


''CMB 
(km s"-') 


(m - M) 
(mags) 


[O/H]^ 
(dex) 


(m - M)"" 
(mags) 


m 
mags 


log 


W 


Inc 
deg 


Type 


Cluster 


IN 


\JZO 


327 


29.80±0.04 


8.55±0.15 


29.81 


8.83± 0.10 


2.328± 


0.024 


65.6 


Sc 




IN Kjyy 




1542 


31.18±0.05 


8.96±0.20 


31.29 


8.32± 0.06 


2.761± 


0.088'^ 


39.0 


Sc 


Fornax 






1400 


31.60±0.05 


9.00±0.15 


31.72 


9.61± 0.07 


2.526± 


0.021=^ 


62.0 


Sc 


Fornax 


NGC 


1637 


671 


30.23±0.07 


9.08±0.15 


30.37 


9.69± 0.06 


2.481± 


0.021 


34.9 


Sc 




NGC 


2090 


994 


30.29±0.04 


8.80±0.15 


30.36 


9.36± 0.10 


2.471± 


0.005 


66.1 


Sb 




NGC 


2541 


696 


30.25±0.05 


8.50±0.15 


30.25 


10.63± 0.14 


2.324± 


0.014 


60.9 


Sc 




NGC 


3198 


877 


30.68±0.08 


8.60±0.15 


30.70 


9.15± 0.10 


2.491± 


0.009 


71.2 


Sc 




NGC 


3319 


979 


30.64±0.09 


8.38±0.15 


30.61 


10.34± 0.13 


2.318± 


0.016 


67.8 


Sc 




NGC 


3351 


1125 


29.85±0.09 


9.24±0.20 


30.03 


8.36± 0.07 


2.568± 


0.036 


43.8 


Sb 




NGC 


3368 


1235 


29.97±0.06 


9.20±0.20 


30.14 


7.82± 0.12 


2.626± 


0.023 


48.4 


Sab 




NGC 


4321 


1897 


30.78±0.07 


9.13±0.20 


30.93 


8.17± 0.06 


2.553± 


0.079 


42.4 


Sc 


Virgo 


NGC 


4414 


981 


31.10±0.05 


9.20±0.15 


31.27 


8.85± 0.05 


2.641± 


0.038 


50.9 


Sc 




NGC 


4535 


2295 


30.85±0.05 


8.85±0.15 


30.93 


9.01± 0.11 


2.585± 


0.037 


45.0 


Sc 


Virgo 


NGC 


4548 


810 


30.88±0.05 


9.34±0.15 


31.08 


8.61± 0.11 


2.634± 


0.040 


35.0 


Sb 


Virgo 


NGC 


4603 


2877 


32.67±0.11 


8.90±0.20'^ 


32.77 


9.85± 0.10 


2.646± 


0.007 


51.0 


Sc 


Gen 30 


NGC 


4725 


1489 


30.38±0.06 


8.92±0.15 


30.48 


7.77± 0.07 


2.625± 


0.030 


62.7 


Sb 


Virgo 


NGC 


7331 


491 


30.81±0.09 


8.67±0.15 


30.85 


7.61± 0.08 


2.732± 


0.006 


64.9 


Sb 





^ 12+log(0/H)'' Corrected for metallicity effects using 
to Will scale using Equation 2. Both galaxies have flat 



the lSakai et al. I l|2004l ) correction"^ Widths from 
ORCs.'* Estimated metallicity as in lNewman et 



ORCs corrected 
aTl 119991) 



ror comes both a combination of the un certain metalhc - 
ity dependence of Cepheid distances (Sa kai et al. 11200^ . 
WFPC2 caHbration systematics, and uncertainty on the 
distance to the LMC which sets the zeropoint for Cepheid 
distances. We account for the metalhcity dependence of 
the Cepheid distance scale using the correction derive in 
ISakai et al. I (|2004l ). Almost all galaxies in the Cepheid- 
SFI++ sample are more metal rich than the LMC, re- 
sulting in an overestimate of Hq if the dependence is not 
accounted. If the metallicity dependence is not corrected 
for we measure Hq = 77 ± 2 km s~^, implying that the 
maximum systematic error which could be introduced if 
this correction is wrong is ^ ±2 km s~^. The effect of 
the uncertainty on the WFPC2 calibration and the dis- 
tance to the LMC are discussed in detail in ISakai et al. I 
()200Cl() where a total systematic error on the Cepheid ze- 
ropoint of 0.16 mag is quoted. This adds a error of ±6 
km to our estimate of Hq, resulting in a final estimate 
of 7?o = 74±2±6 km s'^. 

Our measurement of Hq is identical t o the determi- 
natio n of Hq = 74 ± 2 km s~^ Mpc~^ (jSanchez et al. I 
120061 ) from a combination of WMAP and 2dFRGS data, 
and the recent measurement of Hq = 74 ± 3 ± 6km 
Mpc~^ which uses a new Cepheid calibrati on fr om NGC 
4258 in combination with SNIa distances (|Macri oTI 
|2006[) . Other measurements of Hq using Cepheid cali- 
brations of the TF relation include the HST Key Project 
value of i?o = 71 ± 4 ± 7km s"^ Mpc'^ (|Sakai et al. I 
I2000D . and the SFI measurement of iJo = 69 ± 2 ± 6 
km s~^ (Giovanelli et al. 1997a; note that neither of these 
determinations accounted for the metallicity dependence 
of the Cepheid zeropoint). We point out that the domi- 
nant source of error in our determination of the value of 
Hq comes not from the TF relation, but from the uncer- 
tainty in the zeropoint calibration of the Cepheid rela- 
tion, and ultimately from the uncertainty of our knowl- 
edge of the distance to the LMC. If this error were not 
present TF could measure Hq to a comparable accuracy 
to the best "precision cosmology" available, and better 
than is possible from just WMAP alone which quotes 
Hq = 73 ± 3km s~^ an d also requires an a ssumption 
that the universe is flat (jSpergel et al. 1 120061 ). In fact. 
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Fig. 12. — Galaxies with measured distances from Cepheid vari- 
ables plotted on the TF relation. Sc galaxies are shown as filled 
circles, earlier spirals as open circles. Absolute magnitudes have 
been corrected for morphological type dependence and the bivari- 
ate fit to the in-|- sample relation is shown as the two solid line 
(adjusted for Hq = 74±2km s"'^ Mpc"'^). The bottom panel shows 
the residuals from this relation. Errors on the residuals are the sum 
in quadrature of the measurement errors on the magnitudes, the 
Cepheid distances (including the metallicity dependence) and the 
error on the widths multiplied by the slope of the template (this 
last error dominates). The dotted line shows the error weighted 
mean residual. 

work on imp r oving the Cepheid calibration is underway. 
iMacri et al. I ()2006D discuss the possibility of improving 
Hq measurements based on the Cepheid scale to ±5% 
in the near future, noting that soon there will be four 
galaxies with geometric distance measures (NGC 4258, 
LMC, M31, M33) which could be used to set the Cepheid 
zeropoint. 

8. CONCLUSIONS 

In this paper we present the SFI-I — h sample of spi- 
ral galaxies. Galaxies in this sample have both I-band 
photometry and velocity widths (either from HI global 
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profiles or optical rotation curves) suitable for use in the 
Tully-Fisher relation. From a subset of this sample con- 
sisting of 807 galaxies in the vicinity of 31 nearby clus- 
ters we rederive the I-band TF relation. This sample 
constitute by far the largest ever available to calibrate 
the TF relation. We argue that it provides evidence for 
a type-dependent TF slope that is steeper for later type 
spirals, in effect fitting for two TF relations - one for 
the 465 Sc galaxies in the sample, the other for the 342 
Sa/Sb galaxies. We measure the cluster peculiar velocity 
dispersion for the 31 cluster samples. This quantity is di- 
rectly related to the shape of the initial power spectrum 
of density fluctuations and is used to estimate a value of 
n^ '^as = 0.52 ± 0.06. By comparing the TF template 
with the Cepheid distances to 17 galaxies we measure 
Hq = 74 ± 2 ± 6 km s~^ Mpc~^; an accuracy compara- 
ble to the best precision cosmology. A summary of the 
numeric results is presented in Tabled A more detailed 
summary of the conclusions follows in this section. 

In Section 3.1 we introduced a subset of the SFI-I--I- 
sample which consists of 807 galaxies in the vicinity of 
31 nearby clusters (within vcmb =10,000 km s~-^). This 
sample is divided into an in sample consisting of 483 
bona fide cluster galaxies and an in+ sample which also 
includes galaxies considered to be peripheral members 
of the cluster. Cluster membership assignments are dis- 
cussed further in Springob et al. (2007). 

In Section [312, we discuss various biases which mod- 
ify the observed TF relation from the intrinsic relation 
within a cluster sample. We derive corrections for these 
biases which include morphological type bias and incom- 
pleteness bias due to the implicit magnitude limit of the 
sample. The magnitude of the incompleteness bias de- 
pends most strongly on the assumed scatter in the TF 
relation. 

We find evidence to suggest that the slope of the TF re- 
lation gets shallower for early type spirals and construct a 
width-dependent morphological correction to reconstruct 
the relation for Sc galaxies. This is equivalent to fitting 
two separate TF relations to the sample - one to the ma- 
jor subset of 465 Sc galaxies, a second to the 342 galaxies 
in the sample with types earlier than Sc. The slope dif- 
ference can be explained either by differing completeness 
characteristics for the different morphological types in 
the sample, or a real physical difference in the galaxies. 
We argue for the latter as the difference is shown to be 
largest at the high width end of the relation where the 
bi as correction s are small est. In the numerical models 
of iMo fc Maol (|2000[ ) and iNavarro fc Steinmetel ()2000f ) 
this appears consistent with the idea that earlier type 
spirals have a smaller fraction of their mass in the disk 
at a given rotational velocity, and implies that as well 
as being dimmer than later types at a given rotational 
velocity they have less concentrated halos. 

Any TF relation fit to a sample with an explicit or 
implicit magnitude limit will suffer from biases relating 
to incompleteness. The largest such bias is caused by 
it being possible to observe galaxies which are scattered 
brighter than the TF relation from near the magnitude 
limit, while those scattered dimmer will not make it into 
the sample. We calculate and correct for this bias, as 
well as those relating to the finite depth of the clusters. 
All bias corrections applied are summarized in Section 
3.6. 



Individual TF relations are fit to each of the 31 clus- 
ters samples and presented in Section [H We search for 
environmental dependence of the TF relation by look- 
ing for correlations in global cluster properties with the 
slope, zeropoint and scatter. No correlations are found. 
There is also no dependence of TF residuals on projected 
cluster center difference for galaxies in the sample. We 
therefore argue that there is no evidence for environmen- 
tal dependence of the TF relation. It is noted however 
that because of the morphology-density relation we are 
already implicitly correcting for an environmental depen- 
dence by applying morphological type corrections. 

The individual cluster samples are combined in Section 
Oto create a global TF template. In order to do this, pe- 
culiar velocities are derived for the clusters and a cluster 
velocity dispersion of cr = 298 ± 34 km s~^ is measured. 
This measurement provides information on cosmological 
parameters as it depends on the initial spectrum of den- 
sity fluctuations. The value is completely independent, 
yet within 2a of that predicted from the best fit cos- 
mological model to the third year WMAP data release 
which imphes cr = 239 ± 23 km s~^ (|Spergel et al. | [2006). 
With minimal assumptions about cosmology we use the 
measurement to estimate JI^'^cts = 0.52 ± 0.06 from TF 
data alone. 

Direct, inverse and bivariate linear fits to the com- 
bined in and in+ samples (and subsets of them) are 
presented. Quadratic fits do not reduce the observed 
scatter and so are not reported. We favor the bivariate 
fit to the in+ sample, which gives a template TF re- 
lation of Af - 51og/i = -20.85 - 7.85(logiy - 2.5), or 
L (X v^-^. The statistical error on the zeropoint is ~ 0.02 
mag, an additional 0.03 mag is estimated to account for 
the possible net motion of the cluster reference frame in 
the CMB. The statistical error on the slope is ^ 0.10 
mag""'^. The TF relation obtained here will be applied 
to the rest of the galaxies in the SFI-|--t- sample in fu- 
ture papers including Masters et al. (2007, in prep.) and 
Springob et al. (2007) to study the local peculiar velocity 
field. 

The scatter in the TF relation is as important as 
the relation itself, both in models of galaxy forma- 
tion and evolution and in deriving reliable distance er- 
rors. We measure the total incompleteness bias cor- 
rected scatter from the template relation to be e = 
0.41 — 0.44(log W — 2.5). Once measurement errors have 
been accounted for this results in an intrinsic TF scatter 
of e = 0.35 — 0.37(log — 2.5). This scatter corresponds 
to distance errors of 10% for the largest width galaxies 
increasing to 26% for the smallest widths. 

As an independent measure of the zeropoint, that does 
not rely on the assumption that a subset of the most dis- 
tant clusters are at rest with the CMB, we use 17 galax- 
ies in the SFI-|— I- that also have published distances from 
Cepheid variables. We find a negligible mean offset in the 
zeropoint of the TF relation for these galaxies relative to 
the "basket of clusters" bivariate fit to the in+ sam- 
ple when Hq = 74 ± 2 km s~^ Mpc^^ is assumed, pro- 
viding a reassuring check that the "basket of clusters" 
method provides a reliable measure of the TF zeropoint. 
Turning this method around to estimate Hq from a com- 
bination of the template and Cepheid distances we find 
iJo = 74± 2 ± 6 km s"^ Mpc"^ The dominant source of 
error on this measurement comes not from TF directly. 
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TABLE 6 

Summary of Numeric Results 



Parameter 



Values 



Template relation (bivariate fit to in+ sample) 

Template scatter 

Cepheid calibrated zeropoint 

Hubble's constant 

Cluster peculiar velocity dispersion 



M - 51og?t = -20.85 ±0.03 - (7.85 ± 0.1)(log W - 2.5) mag 

L = 3.7 X 10iO(«;niax/200 km s-'^f-'^LQ/h'^ 

etotal = 0.41 - 0.44(log W - 2.5) mag 

^intrinsic = 0.35 - 0.37(log W - 2.5) mag 

M = -20.49 ± 0.05 ± 0.06 - (7.85 ± 0.1){log W - 2.5) mag 

Ho = 74 ± 2 ± 6 km s"! Mpc"! 

cr = 298 ± 34 km s"! 

QO-Sg-s = 0.52 ill 0.06 



but rather from the uncertainty on the distance to the 
LMC (via the Cepheid cahbration). 

The TF relation stiU has a place in the era of "pre- 
cision cosmology" to provide independent checks of the 
best fit cosmological model. It has the potential to mea- 
sure some parameters (notably the cluster peculiar veloc- 
ity dispersion and Hq) to an accuracy comparable to the 
best available. In addition the TF relation and its scat- 
ter provides valuable information about the formation of 
disk galaxies. However, the reliability of any conclusion 
drawn from TF depends critically on the availability of 
an unbiased template relation, such as the one provided 
here. 
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